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Abstract 

Novel generation of silicon-based photodetectors are attractive alternatives to the traditional phototubes. They offer 
significant advantages but they present new challenges too. Presence of afterpulses may affect many characteristics 
of the photodetectors. Simple statistical model of afterpulsing is used to evaluate the contribution to the observed 
dark count rates, to examine the contribution to the pulse height resolution and to demonstrate the modification of the 
observed timing properties of the SiPMs. 


1. Introduction 

Silicon photomultipliers (SiPM’s, MPPC’s) are the arrays of avalanche photodiodes operating in a Geiger mode. 
Their high gain, high detection efficiency and compact sizes in conjunction with very high pulse height resolution 
enabling ’photon counting’ make them very attractive candidates for replacement of the traditional photomultipliers, 
PMTs, especially in applications where the total area of photodetectors required is small. 

SiPM’s offer several advantages over the PMT tubes but they also present new challenges related to solid state 
nature of these devices. Thermally generated free charge carriers create signals identical to those caused by photoelec¬ 
trons, hence the dark count rates in typical devices are relatively large, in the range 10 4 cps/mm 1 2 - I0 6 cpslmm 2 F\. 

Large number of carriers are present in the Geiger discharge; some of them may be trapped in metastable traps. 
Their subsequent release may produce additional pulses, afterpulses, which contribute to the observed signal. The 
stochastic nature of this process affects the performance of photodetectors. 

Afterpulsing is a result of a large number of independent ’events’ (electron traversing the junction) hence it can be 
well modelled using standard statistical methods. In principle one can derive any observable quantity or distribution 
by a suitable convolution of Poisson and exponential distributions with the mean values and lifetimes appropriate for a 
detector in question. We have calculated these convolutions using Monte Carlo integration and generating 50 million 
initial afterpulses. 

2. Modeling of Afterpulses; Single Pulse 

Afterpulses can be characterized in terms of the probability that a given pulse will produce a subsequent afterpulse 
and their time constant T a f,. In this note the P a f t will be used to denote the average number of trapped electrons 
(hence the number of afterpulses produced by a single pulse, not counting the next generation of afterpulses). For 
Paft 1 it has the meaning of the actual probability of producing an afterpulse. These parameters depend on the 
operating point of the SiPM. Increase of the bias voltage Vbu, s leads to corresponding increase of the P„f, because of 
the increase of the number of charge carriers in the avalanche Q = C ce u(Vbi as - Vbreakdown) = C ce uV ov , where C ce u 
is the cell capacitance and the Vbreakdown is the breakdown voltage, and the simultaneous increase of the probability 
that a free charge carrier will initiate an avalanche, Pceiger- The afterpulsing time constant T a f, may depend on the 
temperature of the photodetector. 

In general there might be several trapping levels present with different trapping probabilities and corresponding 
time constants. In the following we explore a simple model of one trapping level with a single afterpulsing probability 
and one time constant. 
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Afterpulsing is a stochastic process, the actual number of additional pulses undergoes fluctuations according to 
the Poisson distribution with P a f, being the expectation value. In the note the quantity P a f, will be referred to as 
’afterpulsing probability’. Every afterpulse produces a standard avalanche, hence it may lead to a production of 
additional afterpulses. 

In a realistic detector the situation is somewhat more complicated as the bias voltage after an avalanche drops to 
ybreakdown and is being restored with the time constant t«c characteristic for a given device. Electrons released from 
the traps at timescales shorter than t«c have lower probability of producing an avalanche because of reduced Pceiger 
and the charge of the avalanche will be reduced proportionally to the actual value of overvoltage, V ov . The latter effect 
will lead to a reduced afterpulsing probability P a f, of the subsequent afterpulses regardless of their timing. These 
effects depend on the specific devices and they are ignored in the present study. 

Afterpulsing is simulated assuming a single avalanche creating a Poisson-distributed number of afterpulses dis¬ 
tributed with the exponential decay time. The procedure is iteratively applied to each of the generated afterpulses, 
thus leading to trains of pulses with varying numbers of pulses and their time distributions. 


3. Afterpulsing Pulse Rates; Effective Gain 

Afterpulsing with fixed probability would lead to an increase of the number of pulses by an additional contribution 
of 

N aft = Paf , + P 2 afl + P\ ft + ... = (1) 

and it may be treated as an additional gain factor, albeit dependent on the bias voltage. As shown in FigjTjthe same is 
true even in the case when the afterpulses are created with fluctuations governed by the Poisson statistics. 



Figure 1: Average number of afterpulses as a func¬ 
tion of the afterpulsing probability with the actual 
number of afterpulses undergoing statistical fluctu¬ 
ations. 



Figure 2: Excess noise factor as a function of the af¬ 
terpulsing probability. For comparison a linear func¬ 
tion 1 + P a f, is shown as a blue line. 


The additional gain factor undergoes additional fluctuations thus it reduces the intrinsic resolution of the photode¬ 
tector, which can be expressed in terms of the Excess Noise Factor, ENF 1 . Fig. [2] shows that the ENF, calculated 
from the mean value and the RMS of the distribution of the number of afterpulses, grows very rapidly, much faster 
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than linearly, with the P a ft- 

The increase of the number of additional pulses with P a f t (thus with Vbias) is the primary reason for the observed 
increase of the dark count rates with the bias voltage. 

4. Afterpulses Time Distribution 

Afterpulsing time constant can be determined experimentally by studying the observed time distribution of pulses 
following some ’trigger’ pulse. An external light pulse can be used to provide the initial trigger and the contribution 
of thermal pulses can be subtracted. The shape of resulting time distribution is modified by the statistical nature of 
afterpulsing, however. 



time, in units of tau 

Figure 3: Time distribution, in units of trc- of 
the observed pulses for afterpulses probability 
ranging from 0.1 to 0.5. The black distribution 
corresponds to the underlying exponential and 
no subsequent afterpulsing. 



time in units of tau 


Figure 4: Time distribution, in units of trc, of 
the observed pulses for afterpulses probability 
ranging from 0.6 to 0.95. The black distribu¬ 
tion corresponds to the underlying exponential 
and no subsequent afterpulsing. 


Figs. [3 and [4] illustrate the impact of subsequent afterpulsing on the time distribution of the observed pulses 
for different values of the P a f t . The time distribution of the observed pulses retains its exponential shape, but the 
characteristic time constant t 0 /„ grows with the afterpulsing probability, as shown in Fig. [3 

The value of r„/„ is a relevant parameter to describe the performance of a photodetector at the given operating 
conditions but it does not, in general, reflect the time constant of the underlying physics process, r a f t . The latter can 
be determined better by studying the time distribution of the ’first’ afterpulse. Such a method would yield a correct 
result in the case when the number of afterpulses to a given pulse is one. A Poisson nature of the number of produced 
afterpulses leads, however, to a systematic underestimate of the r fl / r as shown in Fig. [6] 

5. Modeling of Afterpulses; Multiple Thermal Pulses 

’Dark’ pulses, i.e. pulses in the absence of the light signal consist of several components. The primary source of 
these pulses are the thermal excitations of charge carriers leading to the avalanche process with the Pceiger probability. 


'Excess noise factor is defined as ENF = —. For purely poissonian fluctuations ENF = 1.0. Any deviation of ENF from unity indicates 
additional fluctuations beyond purely statistical ones. 
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The rate of these pulses depends on the temperature and the bias voltage, but the resulting pulses are randomly 
distributed in time with the overall rate 


R(T, V ov ) = Rthermal(T) X Pceiger(Yo\) (2) 

where R,hermai is a rate of thermally generated carriers at the temperature T. 

The afterpulses (including the subsequent after-afterpulses) will increase the average rate of the observed pulses 
by a factor 1 + N a / r but the time structure of the resulting ensemble of pulses will be depend on the interplay between 
the afterpulsing time constant and the typical time distance between the thermal excitations A t t i ierma i. The dark pulses 
will be distributed randomly in time when A t,i ierma i <sc T a f, or there will be randomly distributed trains of pulses (with 
the train multiplicity and duration dependent on P a ft) when At thermal r a / f . 

Thermally produced pulses have been generated in the interval 0 - 10 /jsec and the resulting afterpulses have been 
generated with the average multiplicity given by P a f, and time constant r a / t = 50 nsec. Fig. |7]shows an example of the 
time distributions of pulses simulated with the thermal random rate of 2 MHz and different afterpulsing probabilities. 
The apparent raise of the instantaneous rate is a reflection of the lack of afterpulses corresponding to thermal pulses 
produced at t < 0. the time interval required to reach the asymptotic level in such a situation depends on the interplay 
of the afterpulsing time constant and the afterpulsing probability. 

6. Thermal Pulses Rate Estimates 

In the presence of afterpulses the observed ’dark’ pulse rate is an overestimate of the rate of random thermal 
excitations by a factor 1 + N a f t which can exceed an order of magnitude. A better estimate of the initial thermal rate 
can be obtained from a distribution of pulses multiplicity in some time interval At. The probability of observation of 
zero pulses P( 0) in the given time interval provides an estimate of the average number of pulses expected in this time 
interval 

< N >- -logP( 0) (3) 

If there were no thermal pulses in the specific time interval then there will be no afterpulses either, therefore such 
method appears to be systematically superior to the use of the raw observed rate. A possible presence of the afterpulse 
to the pulses preceding the time interval in question reduces the probability of observing zero pulses and biases 
the resulting < N > towards the higher values, however. To study the resulting bias the random thermal pulses 
corresponding to the rate of 2xl0 6 cps were generated in the time interval of 10 jisec with afterpulses corresponding to 
the time constant r a / f of 50 nsec. To avoid an additional bias due to the imperfection of the simulation at the beginning 
of the simulated time interval the method was applied using variable gate length, all of the gates starting at t = 7.5 /usee, 
well inside the steady state even for large afterpulsing probability, see Fig. [7] 

Figs. [8] and [9] illustrate the magnitude of the systematic bias introduced by afterpulsing in the determination of 
the thermal pulses rates as a function of the of used gate length and the afterpulsing probability. For the gate length 
longer than 1 psec — 20 r a ft the true thermal rate can be determined with the accuracy better than 20% even for very 
large afterpulsing probability. 


7. Afterpulsing and Absolute Calibration 


Absolute calibration of photodetectors, i.e. conversion of the observed signals into the number of photons impact¬ 
ing the photodetector, is the crucial step in several classes of experiments. One of they methods successfully used with 
photomultipliers involves illumination of the photodetector with a light signals of constant intensity and determination 
of the average number of photons impacting the photodetector from the width of the observed pulse height distribution 
using the Eq. 0] 


< N >= 


Mean 2 
RMS 2 


(4) 


where Mean is the mean value of the measured distribution of signals and RMS is the width of this distribution. This 
method relies on the assumption that the variation of the observed signal is dominated by the Poisson fluctuations of 
the number of detected photons. It is a very good approximation for the PMT tubes, which have ENF very close 
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to one. Fluctuations of afterpulsing contribute significantly (Fig. [2} to the fluctuations of the observed signals thus 
reducing the ’determined’ number of the incident photons with respect to the true size of the photon signal. This effect 
can be parameterized as 

< N det >=< N ,rue > (1 - P aft ) (5) 

where < N de! > is the number of incident photons determined using this method whereas < N ,rue > denotes the 
number of photons and it is illustrated in Fig. [TO] 


8. Summary 

Afterpulsing is an intrinsic feature of SiPM-type photodetectors and it does contribute to the various factors rele¬ 
vant to the performance of these photodetectors. The afterpulsing probability is proportional to the overvoltage, hence 
many of the observed characteristics of the SiPMs depend on the bias voltage. 

A simplified model of afterpulsing gives and explanation of a rapid increase of dark count rate (or a dark current) 
with the bias voltage and, at the same time, it predicts a serious deterioration of the intrinsic resolution of the detectors. 
The observable timing characteristics of SiPMs depend on the bias voltage as longer and longer trains of afterpulses 
are generated with the increasing voltage. Various methods of the experimental determination of the true underlying 
dark count rates and/or absolute calibration are affected by the afterpulsing and therefore their outcome depend on the 
operating point of the devices. Knowledege of afterpulsing level can be used to derive the corrections to determine 
the underlying values of these parameters. 

Proper inclusion of the afterpulsing in the analysis of the data is very important, but it depends on the properties 
of the photodetector (including the overall rate of afterpulsing and its time constant as well as the RC recharge time 
of the cells) as well as on the details of the experimental analysis (gate length). In the simplified model described 
above the afterpulsing rate is identical with the rate of thermally induced excitations. In the realistic detector the 
rate of afterpulses will be reduced by the combine effect of Geiger probability and the recharge of the cell, therefore 
the relevant P a f, will be in general much smaller then the thermal excitation rate and the effective time constant 
of afterpulses will be diffrent from the actual lifetime of trapping centers. Nevetrheless it is expected that using the 
effective P a f, and the corresponding time constant should provide a good description of the real devices and the results 
presented here can be used as a guidance to derive the corrections relevant for specific applications. 
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Figure 5: Effective afterpulsing time constant, in 
units of trc, as a function of the afterpulsing prob¬ 
ability. 
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Figure 6: Time constant, in units of trc, deter¬ 
mined from the ’first afterpulse’ as a function of 
the afterpulsing probability 
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Figure 7: Instantaneous rate as function of time after a ’turn-on’ for different afterpulsing probabilities. The thermal 
random rate is 2x 1 (f'cps and afterpulsing time constant is 50 nsec. 
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Figure 8: Ratio of the rate of thermal pulses es¬ 
timated from the P( 0) to the true rate of the un¬ 
derlying thermal pulses for different afterpuls¬ 
ing probability P a f t as a function of the length 
of the time interval used. 



Figure 9: Ratio of the rate of thermal pulses 
estimated from the P{ 0) to the true rate of the 
underlying thermal pulses for different time in¬ 
tervals as a function of afterpulsing probability. 



Figure 10: Reduction of size of the photon signal using the Eg|4]with respect to the truth as a function of Afterpulsing 
probability. 
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